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Stellar model fits and inversions 
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The recent asteroseismic data from the CoRoT and Kepler missions have provided an entirely new basis for investigating 
stellar properties. This has led to a rapid development in techniques for analysing such data, although it is probably fair to 
say that we are still far from having the tools required for the full use of the potential of the observations. Here I provide a 
brief overview of some of the issues related to the interpretation of asteroseismic data. 
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1 Introduction 

The goals of the asteroseismic analysis are evidently to im- 
prove our knowledge about stellar properties and our un- 
derstanding of stellar structure and evolution. More specif- 
ically, simple model fits, or even direct use of scaling 
relations, provide estimates of the masses and radii of 
stars, while more detailed analyses, including also other 
observations of the stars, yield more precise values of 
mass and radius as well as an estimate of the stellar age. 
Such determinations clearly have substantial broad inter- 
est, not least for stars known to be hosts of planetary 
systems. However, from the point of view of stellar as- 
trophysics we clearly wish more information, and indeed 
we can go much further from the data now becoming 
available, not least from the CoRoT and Kepler missions 
dMichel et al.ll2008l:jBaglin et alJl2009l:lBorucki et alj |20iri 
Gilliland et al.l2010tlKoch et al.l2010l) . Recent or upcoming 
highlights are the inference s on the core structure and rota- 
tion of red-giant stars (e.g.JBedding et alj|201~ ; Beck et al 



2012; Mo sser et al.ll2012 a b) and the detection of effects in 
solar- like oscillations of acou stic glitches in main-sequence 
stars ( Mazumdar et al J[20 1 2.allbb . Fortunately, both missions 
have now been extended until 2016. 

Much of the data fitting or other form of analysis is done 
within the framework of stellar models characterized by a 
specific set of parameters. A crucial goal of asteroseismol- 
ogy is to test whether such models adequately represent real 
stars and, if not, how the model physics should be improved. 
Inadequate modelling should be reflected in a failure of the 
models to fit the observed stellar properties, particularly 
the frequencies, to within their observational errors, for any 
choice of parameters characterizing the model. At this point 
one then has to determine how the models should be further 



improved, perhaps through inverse analyses to try to local- 
ize the cause of the discrepancy. It is evident that the identi- 
fication of such significant departures requires great care in 
the statistical analysis of the data, and it is probably fair to 
say that we have so far not quite reached this point. 

In this brief review I concentrate on the determination 
of stellar structure based on solar-like oscillations, with 
main emphasis on stars on or near the main sequence. 
However, it should be noted in passing that striking re- 
sults have also been obtained based on CoRoT and Ke- 
pler data fo r more massive stars sh owing heat-engine driven 
pulsations (Degroote et al. 2010) and for the subdwarf B 
stars (e.g. JVan Grootel et al .2010). Also, extensive ground- 
based campaig ns have yielded fund amental results for white 
dwarfs ("e.g.. iMetcalfe et alj|2004|) . Thus asteroseismic in- 
vestigations are covering a very broad range of stellar prop- 
erties and stellar evolutionary stages. 

2 Properties of stellar oscillations 

As a background for the discussion, I present a brief 
overview of the relevant properties of stellar oscillations. 
Much more detailed descript ions were provided, e .g., by 
Christensen-Dalsgaard! d2004 and lAerts etal] d2010b . 



2.1 Basic properties 

We consider only oscillations that are solar-like, in the 
sense that they are intrinsically damped and excited by 
the acoustic noise generated by near -surface convection, 
which reaches near-sonic sp eeds (e.g. jGoldreich & Keelev 



1977t iHoudek et al.l 1 19991) . This mechanism depends 
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on the near-surface properties of the oscillations and 
the time-dependence of t he most vigorous convection 
( Kieldsen & Beddfngll201 lb . making it most efficient at fre- 
quencies approaching the acoustic cut-off frequency in the 
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stellar atmosphere, approximately given, in terms of cyclic 
frequency, by 

GM 



Cphot 



D2T 1 / 2 
- fl J cff 



(1) 



where c p hot and H p are the sound speed and pressure scale 
height at the stellar photosphere. The second expression as- 
sumes the ideal gas law and that the atmospheric tempera- 
ture scales like the effective temperature Toff; here G is the 
gravitational constant, and M and R are the mass and radius 
of the star. 

In relatively unevolved stars, on the main sequence, the 
modes in the frequency range of efficient driving are acous- 
tic modes, with cyclic frequencies that are asymptotically 
given by 



v n i ~ Av{n + 1/2 + e)-5 



(2) 




dVandakurov et alj Il967t iTassoull Il98d: iGough] 1 1986b . 
where n and I are the radial order and spherical harmonic 
degree of the mode. Here Av is approximately the inverse 
acoustic diameter of the star, 



(3) 



where c is the adiabatic sound speed and r is the distance 
to the centre of the star. As expected from homology argu- 
ments this scales approximately as the square root of the 
mean density of the star, 

/M\ 1/2 

Also, e is a phase depending on the near-surface proper- 
ties of the mode and 5qi is a small correction that in main- 
sequence stars is sensitive to the central regions of the star. 

In terms of the observed frequencies the asymptotic 
properties in Eq. (f2j) are often utilized by considering the 
large frequency separation 



Av n 



-1 / 



and the small frequency separations 

Sv n i = V n l — Vn-ll+2 ■ 



(6) 



Here Av n i approximates Av, although with some varia- 
tions with frequency that largely arise from the near-surface 
acoustic glitches (see Section [3~3l : 5v n i, as Sou provides a 
measure of conditions in the core of the star and hence, for 
main-sequence stars, the stellar age. 

The amplitudes of the modes are determined by a bal- 
ance between the energy input from the stochastic excitation 
and the damping. This leads t o a bell-shaped distrib ution 
of amplitudes with frequency dGoldreich et al J 1 19941) with 
a typically relatively well-defined maximum at a frequency 
v mM . There is observational evidence that v max scales as 



v ac , i.e., 



GM 



R 2 T 



1/2 



(7) 



(e.g.. iBrown et all Il99ll; iBedding & Kieldsenl 12003 



Stello et al.ll2008l) . with some although perhaps not full the- 
oretical justific ation in terms of the near-surface properties 
of the modes (Belk acem et al ■11201 lh . 

The large and small frequency separations and ^ max 
form the basis for ensemble studies of large samples of stars 
and hence are typically determined by automated pipeline 
analysis method s. A careful comparis on of such techniques 
was provided bv lVerner et al.1 d201 lab . 

As the star evolves beyond the main sequence, the core 
contracts strongly and hence the internal gravitational accel- 
eration becomes very large. Furthermore, hydrogen burn- 
ing leaves behind strong gradients in the chemical composi- 
tion. These effects lead to a marked increase in the so-called 
buoyancy frequency N, which for an ideal gas can be ap- 
proximated by 



N 



g2p (X7 
(V a d 

P 



v + V, 



(8) 



where p is pressure, p density, g is the local gravitational 
acceleration, V = dlnT/dlnp, T being temperature, V a( j 
is the adiabatic value of V and V M = dln/z/dlnp, where 
[i is the mean molecular weight. N is the characteristic fre- 
quency for internal gravity waves, and when N is large the 
star has mixed modes, which behave like acoustic modes in 
the outer parts of the star and as gravity modes in the deep 
interior. For sufficiently high N, found in sub-giant and red- 
giant stars, such mixed modes have frequencies as high as, 
or exceeding, v ac and hence are excit ed efficiently by th e 
stochastic near-surface processes (e.g. JPupret et alj[2009l) . 
Given the dependence of N on the structure of the deep in- 
terior of the star, including the composition profile, these 
modes are powerful diagnostics of the stellar interior. 

High-order asymptotic g modes are uniformly spaced in 
period n n /, which asymptotically satisfies 

n 



(n + e g ) , 



(5) where 



n = 2tt 2 



N- 



dr 



(9) 



(10) 



the integral being over the region where the buoyancy fre- 
quency exceeds the oscillation frequency, and e g is a phase. 
For red giants this leads to a very dense spectrum of high- 
order g modes interacting with the acoustic modes in the 
region of stochastic excitation. For modes with I = 1 sev- 
eral such mixed modes are typically observed in the vicinity 
of each acoustic mode. This allows a determination of the 
corresponding period spacing Alii = Hq/ a/2 which, ac- 
cording to Eq. ( [Tol l, i s a diagnostic of the deep stellar inte- 
rior where N is larg e dBeck et al.l201 1 ; Bedding et al .2011 
Mosser etai]|20 111) . 



off 



In a spherically symmetric star the frequencies do not 
depend on the azimuthal order m of the spherical harmonic. 
Rotation lifts this degeneracy, giving rise to a rotational 
splitting which depends on the internal rotation of the star, 
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weighted by the eigenfunction of the mode. Observations 
of the rotational splitting in the Sun have provided detailed 
inf ormation about the sola r internal rotation (for a review, 
see iThompson et al.ll2003l) . and Kepler has provided strik- 



ing observations o f the splitting of solar-like oscillations in 
evolved stars (e.g., iBeck et al .2012; Dehe uvels et al 

EH! 



Mosser et al. 201 2bb . However, here I ignore rotation and 
concentrate on the determination of stellar properties based 
on the average multiplet frequencies v n i . 

2.2 Near-surface problems 

Modelling of near-surface layers of a solar-like star is com- 
plicated by the presence of convection. The temperature 
structure of the model is determined by the treatment of 
convective e nergy transport, often using the mixing-length 
formulation (IBohm-Vitensel 1 1 9580 . and the models nor- 
mally ignore the dynamical effects of convection, in the so- 
called turbulent pressure. Also, the frequency calculations 
are most often carried out in the adiabatic approximation, 
even though the oscillations are strongly nonadiabatic in the 
superficial layers. However, the treatment of the nonadia- 
batic effects needs a description of the interaction between 
convection and pulsations, through the perturbation to the 
convective flux; also, the perturbation to the turbulent pres- 
sure may affect both the frequencies and the stability of the 
modes. An ove rview of these complications was given by 
Houdekl d2010h . It should be noticed that the uncertain as- 
pects of the structure and oscillation modelling are all con- 
centrated very near the surface. Thus they have little effect 
on modes of low frequency which are evanescent in this re- 
gion, with a very low amplitude compared with the ampli- 
tude in the bulk of the star. 

These near-surface problems dominate the differences 
between observed so lar oscillation frequencies and frequ en- 
cies of solar models ( Christensen-D alsgaar d et al.lll996 ). In 
this case it is possible to separate the near-surface effects 
from the differences between the Sun and the model in the 
deeper parts of the star. A convenient way to do so is to use 
the differential for m of the Duvall asymptotic expre ssion 
for the frequencies (iChristensen-Dalsgaard et alJI 19891) . ac- 
cording to which the frequency differences 5v between the 
Sun and a model satisfy 

s nl — ^n 1 ^ n i/L) + n 2 {iy n i) , (ii) 

V-al 

where L = + 1) and S n i is a scale factor which can be 
calculated from the model, and which may be chosen to be 
close to 1 for low-degree modes. The term in Hi depends 
on the location of the inner turning point of the mode, de- 
termined by v/L, and reflects the difference in sound speed 
between the Sun and the model throughout the star, whereas 
% contains the contribution from the near-surface region. 
In accordance with the above comment the arbitrary addi- 
tive constant in the definition of H2 may be chosen such 
that H-2 is zero at low frequency. In Fig. Q] the solid curve 
shows v'K'xiy) determined in this manner, corresponding to 



the frequency shift 8v caused by the surface effects in the 
Sun. 

In principle, these effects must be taken into account in 
any analysis of frequencies of solar-like oscillations. They 
affect not only the individual frequencies but also frequency 
combinations such as the large and small frequen c y sepa- 
rations (cf. Eqs [5] and [6}. iRoxburgh & Vorontsovl (I2003al) 
pointed out that separation ratios, such as 

f'02 = , (12) 

Vnl — V n -l 1 

can be constructed which are es sentially independent of 
the near-surface problems (see also Otf Flo ranes et al.l2005 



Roxburgh! 120051) . These combinations therefore provide 



clean diagnostics of the core properties of the stars. 

To fit individual frequencies it is necessary to correct 
for the near-surface effects. For distant stars, where only 
low-degree modes are observed, an unconstrained determi- 
nation of the f r equen cy corrections cannot be carried out. 
Kiel dsen et al.l d2008l) suggested to base the correction on 
the solar results, fitting the observed frequencies to 



Vnl 



(best) 

y , 



u (best) ( 

where ui, ' = rvi 



Vnl 

f) 



(13) 



, are the frequencies of the best- 
fitting model, taking into account the near-surface correc- 
tion, and v^ fS> are frequencies of a computed reference 
model, assumed to be close to the best model. The scaling 
with the constant factor r assumes that the frequencies of the 
models are related by homology. In Eq. dT3l the exponent b 
is obtained from fits to solar frequencies, and vq is typically 
taken to be an average of the observed frequencies, which is 
generally close to f max . Kjeldsen et al. provide a procedure 
for determining the amplitude a of the correction and the 
frequency scale factor, by fitting the mean large frequency 
separation and average frequency for radial modes. 

Math ur et alj d2012) carried out detailed fits to a num- 



ber of stars observed by Kepler, including the Kjeldsen et al 



(2008) correction. They found that the amplitudes a roughly 
scaled as Av. According to Eq. (0 this corresponds to a 
roughly consta nt average shift in t he phase e. A similar shift 
was found by White et al. ( 201 lab in a comparison of fitted 
values of e from observed and model frequencies. This sug- 
gests that a more detailed analysis of the phase, including 
how it is affected by the near-surface layers, may provide a 
deeper insight into the properties of the frequency correc- 
tion. 

Although the power-law correction assumed in Eq. dT3l > 
has some theoretical supp ort in an analysis of simple near- 
surfac e modifications ( Christensen-D alsgaar d & Goughl 
1980) it is clearly only an approximation to the frequency 



differences obtained for the Sun (see Fig.Q]). An alternative 
is to base the correction more closely on the functional form 
obtained in the solar case. It is perhaps not unreasonable 
to assume that the relevant frequency scale is the acoustic 
cut-off frequency v ac and hence to replace Eq. dT3l by 

v m = ^° St) + aQ Q {v n i/v ac ) , (14) 
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Fig. 1 In panel a) the solid line shows the surface con- 
tribution to the differences between the solar and model 
frequencies, inferred from analysis of modes over a broad 
range in degree. The symbols show the result of fi tting this 
function to BiSON frequencies (see lChaplin et al.ll2.Q02h of 
degree I = 0—3. For comparison, the dashed line is the cor- 
responding fit of a power law (cf. Eq. [T3j with exponent 
b = 4.9 In panel b) the symbols show the the differences 



betwe en observed frequencies of 16 Cyg A (Metcalfe et al 



2012) and frequencies of a fitted model, without the surface 
correction that was used in the fit. That surface correction 
is shown by the dashed line, while the solid line shows the 
scaled solar surface contribution from panel a), assumed to 
be a function of the frequency in units of the acoustic cut-off 
frequency. 



where the function Q@ is determined from vHiiy), illus- 
trated in Fig. Q] as obtained from fitting solar frequencies 
over a broad range of degrees. The scale factor a and the 
frequency scale r can be determined from a least-squares fit 
to the differences between the observed and the reference- 
model frequencies. In Fig.[T]the symbols illustrate the result 
of applying this procedure to observed solar frequencies of 
degree / = 0—3, corresponding to what can be observed 
in distant stars, and using the same reference model as was 
used to determine Qq . Here the scale factor was extremely 
close to 1 . For com parison, the dashed cu rve shows the re- 
sult of applying the IKieldsen et al.l (12008b procedure to the 



same set of frequencies, with the same reference model. It 
is evident that, while the power law reproduces the general 
rapid increase of the effect with frequency, the detailed fit to 
the extensive and very accurate solar frequencies is some- 
what questionable. 

As a more realistic illustration of the fit to the sur- 
face effect Fig. [TJ5 shows results fo r 16 Cyg A, based o n 
the Kepler observations presented bv lMetcalfe et al.l (|2012). 
I made a fit of the observed freq uencies to stellar mod- 
els, using the lKieldsen et al.l (|2008) surface-term procedure, 
and subsequently determined the best-fitting surface cor- 
rection according to Eq. (fT4l . The dashed and solid lines 
show the power-law surface term and the term correspond- 
ing to Eq. (fT4l . respectively, while the symbols show the 
differences between the stellar and (uncorrected) model fre- 
quencies. Here the advantage of using the solar fit is less 
dramatic, but certainly still significant. Given that 16 Cyg 
A is quite similar to the Sun (with M = l.lM© and 
T c g = 5825 K) that is probably not surprising. However, 
it is interesting that data are now becoming available of suf- 
ficient quality to test more detailed properties of the surface 
effects in stars other than the Sun. Also, the results suggest 
that there might be an advantage in using the solar surface 
functional form in fits to observed frequencies, rather than 
the power law, at least for stars in the vicinity of the Sun. In 
any case a comparison of the results of using these two rep- 
resentations might provide some measure of the systematic 
error introduced in the fit by the use of the surface correc- 
tion. 

It is obvious that the assumption underlying these near- 
surface corrections, namely that they resemble the solar 
correction, is in general questionable, increasingly so for 
stars with propertie s further from those of the Sun. In fact 
Dogan et al.1 ( 120 101) were unable to find a solar-like cor- 



rection for the F5 star Procyon A. Thus some care is re- 
quired in applying the correction, and further tests are re- 
quired to investigate whether it might lead to systematic 
errors in the outcome of the fits to the observed frequen- 
cies. It is evident that the broad range of stars observed 
with Kepler provides an excellent basis for such tests. In 
the longer run the goal must obviously be to improve the 
modelling of the near-surface layers. In the solar case some 
improvement has resulted from the use of detailed hydro- 
dynamical models of the stru cture of the outermost lay- 
ers (e.g.. iRosenthal et al. I l 1999b . while the success of using 
time-dependent convection modelling in the frequ ency cal- 
culations has so far been som ewhat questionable (IHoudek 
l2010tlGrigahcene et alJExTj . 



3 Fitting stellar models 

3.1 Using scaling relations 

The simplest analysis of asteroseismic data is based on the 
overall properties of the oscillations, typically the large sep- 
aration /S.v and the frequency u max of maximum power. Us- 
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ing the scaling relations, Eqs (01 and @, normalizing to so- 
lar values Az/ Q and v max q, we immediately obtain 



R 

R Q 

M 

Mo 



^max,0 



'max,© 



( T e g 



Av Q 

At/ 
A^o 



V^eff,i 



1/2 



3/2 



off,© 



(16) 



(Kall inger et alJuOlOatlMosser et al.ll2010l) . Assuming that 
T ff has been determined, this provides values of R and M, 
independently of any models. However, it should be noticed 
that Av and ^ max appear to quite high powers, particularly 
in the expression for the mass, and hence the results are sen- 
sitive to errors in the observed quantities. 

It is obvious that such a determination of the stellar mass 
and radius depends on the validity of the assumed scalings. 
Since the homologous relation between the relevant stars 
is not exact, we may expect departures from the scaling of 
Av. This was investigated by White et al. ( 201 lbh . for stars 
on the main sequence and the ascending red-giant branch. 
They found a variation of a few per cent which was mainly 
a function of effective temperature; since a significant de- 
parture from homology is caused by the varying depth of 
the convective envelope, which is largely determined by the 
effective temperatu re, such a dependence is not unexpected. 
Miglio et alj d2.012l) found a difference of about 3 per cent 
between Af in models of ascending red-giant stars and he- 
lium burning clump stars at given mean density, related to 
the differences in their internal structure. It should also be 
kept in mind that the near-surface problems discussed in 
Section 12.21 have an effect on Az/, given their strong fre- 
quency dependence. This may cause systematic errors in the 
results of the fit. 

The scaling of i^ max , Eq. @, has no fully secure the- 
oretical basis, since it is not yet possible to make reli- 
able predictions of the amplitudes of stochastically excited 
modes or their dependence on frequency. However, as men- 
tioned above, the relation has substantial observational sup- 
port. A crucial test of this so-called direct technique for 
analysing the asteroseismic data, or other similar asteroseis- 
mic techniques to determine the stellar radius, is evidently 
comparison to independent determinations from the stellar 
angular diameter and distance. Such te sts appear to con- 
firm the validity of the sc aling relations dHuber et alJl2012 ; 
Silva Aguirre et al ]l2012h .A similar test of mass determina- 
tions for individual stars has apparently not been possible. 
However, application to red-giant stars in open clusters pro- 
vides a tes t that the resulting masses span the expected small 
interval dBasu et alj|201 UlMiglio et al.ll2012l) . 

A more strongly constrained, and hence potentially 
more precise, determination of stellar parameters can be ob- 
tained by fitting the observed asteroseismic quantities, to- 
gether with the effective temperature and composition, to 
grids of stellar models. The determination of Av and v max 
for the models in the grid can be based either on the scaling 
relations or, in the c ase of Av, on computed frequencie s 
for the models (e.g.. Istello et alJl2009t iBasu et all [2010). 



The accuracy of the results is clearly sensitive to the stel- 
lar models, and hence it is important to test the effect on the 
fits of using grids computed with d ifferent evolution codes 



(15) dGai et alJ 1201 lb IBasu et alJ 120121) . The grid-based fitting 



was ext ended to include also the small frequency separa- 
tions by lOuirion et al. (2010) in the so-called SEEK algo- 
rithm, which also included a statistical analysis of the re- 
sults. 

These techniques provide a fast and potentially auto- 
mated way to obtain basic stellar properties from the aster- 
oseismic data, and hence they are well suited to the analysis 
of large samples of stars, in what has been called ensem- 
ble asteroseismologjm This has been extremely valuable 
for stellar populatio n studies for main-seq uence stars (e.g., 



Chaplin i 



ir population studies tor main-sequence stars (e.g., 
et alJ l201lb IVerneret alJ 1201 lbh a n d red giants 



2010j; iHekker etal. 



^Kallinger et alJ I2010bt iMosser et al. 
1201 1 alibi iMiglio et alJl2012l) ~ 



3.2 Fits to individual frequencies 

The basic asteroseismic quantities Av and i/ mlK can be 
determined from data with even rather poor signal-to- 
noise levels and, as discussed above, may provide a de- 
termination of stellar global properties. However, with 
the long timeseries obtained from CoRoT and Kepler it 
is oft en possible to determine individual frequencies v n i 



is onen possiDie 10 ueiermine muiviuuai iiequencies v n i 
(e.g.. lAppourchaux et all 120081 120121: iBarban et al.l l2009t 
Campante et al.ll201 li llvlathur et al.ll201 lh . This obviously 
provides far more information about the stellar properties. 

The analysis is typically carried out by fitting the fre- 
quencies to a set of models, characterized by parameters 
{'Pk} which typically include the mass, age and initial 
chemical composition of the star as well as, possibly, the 
mixing length and a parameter characterizing convective- 
core overshoot. In its simplest form the analysis aims at 
minimizing 



(obs) 



(mod) v 



xl = 



N — 1 



cr{v n i) 



(17) 



as a function of the parameters {Vk}- Here N is the num- 
ber of observed frequencies v^ s \ v^' 



are the model 



frequencies, possibly corrected for near-surface effects ac- 
cording to Eq. ( fT3l and <j{v n i) are the standard errors in the 
observed frequencies. To this may be added a similar term 
from other observed properties of the star, such as effective 
temperature, surface gravity, [Fe/H], and, when available, 
radius and luminosity, to obtain the \ 2 which is minimized. 

As is common in multi-parameter fits, \ 2 often shows 
multiple local minima, and a blind search risks landing 
in what is not the globally smallest value. Genetic algo- 
rithms provide a relatively effici ent technique for local- 
izing the global minimum (e.g., iMetcalfe & Charbonneau 
2003). An implementation for the analysis of solar-like 



1 although, according to Douglas Gough, sinasteroseismology might be 
a more appropriate term 
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oscillations was presented by iMetcalfe et al. ( 20091) : this 
also included a local analysis based on singular-value de- 
composition (SVD) to locate more precisely the true min- 
imum, based on the typically relatively coarsely sam- 
pled genetic fit, a nd determine the statistical properties of 



In evolved stars the observation of mixed modes pro- 
vides much more stringent constraints on the stellar param- 
eters. The buoyancy frequency (cf. Eq. [8]) changes rapidly 
with age and so, therefore, do the frequencies of modes 
with a strong g-mode component. This provi des a very pre 



the fit (see also iBrown et al.l 119941) . This has been im- rise d etermination of the stellar age (e.g., IMetcalfe et al 



plemented in the AMP packag^l for analysis of solar- 
like oscillation frequencies, which has be en made avail 
able to the asteroseismic community dWoitaszek et al. 



20101) . This is based on the Aarhus stellar evolution code 
(ASTEC; Christensen-Dalsgaardl2008a ) and pulsation code 
(ADIPLS- Tchristensen-Dalsgaardboo'ibl) . 



Although the genetic algorithm is effective in localizing 
the absolute minimum in x 2 , there may be parameters which 
provide a fit that differs insignificantly from this best fit. In 
some cases this may correspond to p arameters rather dif - 

EoToh . 



Metcalfe et al 



ferent from those of the best fit (e.g.. 
Such non-uniqueness, which is a clear indication of the in- 
trinsic nonlinearity of the dependence of the frequencies on 
the parameters, requires further study. It is likely that it may 
be broken by increased precision in the observations, an in- 
creased range of the observed frequencies or possibly the 
inclusion of additional observables in the fit. 

An unconstrained fit may yield optimal fits that can be 
judged to unphysic al from independent information. An ex- 
ample is the fit by Mat hur et all d2012l) . using AMP, to 22 
Kepler stars; here the initial helium abundance Yq result- 
ing from the fit in several cases was substantially below the 
cosmologicall y determined p rimordial helium abundance 
Yp = 0.248 dSteigmanll20071) . In some cases a value was 
obtained at the lower limit of the range in Yq allowed in the 
fit. Although exceptionally interesting, if correct, the low 
value of Yq is more likely caused by problems in the fitting 
or the underlying models. A related, if less evident, exam- 
ple could be an inferred Yq and initial heavy-element abun- 
dance Zq which depart strongly from the normally assumed 
relation from Galactic chemical evolution. Also, for an oth- 
erwise sun-like star a value of the inferred mixing length 
very far from the solar calibration might be suspect. 

Potential problems of this nature can be avoided by in- 
cluding prior information in the fit, although one must ob- 
viously be careful not to let the result be too strongly influ- 
enced by perhaps unwarranted prejudice. Such prior infor- 
mation can be incorporated through a Bayes ian approach. 



as w as also done in the SEEK algorithm dOuirion et al 



2010h . An interesting example, applied to fits of individ- 
ual frequencies to grids o f stellar models, was provided by 
Gruberbauer et al.l (12.012b . They demonstrated that this pro- 
cedure can also be used to suppress the effect of the sur- 
face problems. A statistical characterization of inferred stel- 
lar parameters, usi ng also Bayesian priors, was obtained by 
Bazot et alJ (120121) using a Markov Chain Monte Carlo al- 
gorithm in a fit to observations of a Cen A, although at very 
considerable computational expense. 



2010)0 A powerful technique for the analysis of da t a with 
mixed modes was presented bv lDeheuvels & Michell(l201lh 
and applied t o CoRoT data (for a n appl ication to Ke- 
pler d ata, see iDeheuvels et al. I l2012l) . Also, Benomar et al 



d2012l) showed how additional information could be ob- 
tained from the detailed properties of the frequencies of 
such mixed modes. 

A major breakthrough in the seismic analysis of red gi- 
ants has been the detection of mixed dipolar modes with a 
strong g-mode componen t ( Beck et ap201 1 ; Bedding et al 



mg 

201 It iMosser et al.11201 11) . These appear as a 'dipolar for- 



est', a group of peaks around the position of the dipolar 
acoustic resonances where the modes have their strongest 
p-mode character. The period spacings between the modes 
in such groups reflect the uniform period spacing for pure 
g modes (cf. Eq. |5), although with departures caused by 
the interaction with the acoustic cavity. It is possible to 
extrapolate the observed period spacings to determine the 
true g-mode spacing Alii and hence constrain the prop- 
erties of_Ae buoyancy frequency i n the core of the star 
(e.g., IStelloll2012l) . It was found by iBedding et al.l d201ll) 
that AITi provides a clear distinction between stars on 
the ascending red-giant branch, where the energy comes 
solely from the hydrogen-burning shell around an inert he- 
lium core, and the 'clump stars', where in addition the re 



is he lium fusion near the centre (see also IMosser et al 



201 1). This distinction comes in part from the fact that the 
helium-burning stars have a small convective core, where 
consequently the buoyancy f requency is essentially zero 
( Christensen-Dalsgaardll2012a ). An asymptotic expression, 
characterizing the frequencies in terms of the period spac- 
ing and the coupling between the gravity-wave and acoustic 
regions, was derived by Gou pil (in preparation ), based on a 



general asymptotic analysis (Shibahas M|1979[ lUnno et al 



19791) (see also Christensen-Dalsgaard 2012b). This was 
used by Mosser et al.l ( 2012ah for an ensemble analysis of 
red giants observed by Kepler. Procedures for detailed fits to 
the observed frequencies of mixed dipolar modes are under 
development, and will probably benefit from an application 
of this asymptotic expression. 

3.3 Analysis of acoustic glitches 

In the analysis of stellar oscillation frequencies informa- 
tion about specific features of the star may be obtained 
through suitable combinations of frequencies, often based 
on asymptotic analyses of the properties of the oscilla- 
tions. An important example is the effect of sharp features, 



Asteroseismic Modeling Portal; see http : / / amp . ucar . edu/| 



It should perhaps be emphasized that the accuracy of such determina- 
tions still depends sensitively on uncertainties in the stellar modelling. 
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Fig. 2 Effects of acoustic glitches on observed solar fre- 
quencies. In the top panel the symbols show second differ- 
ences obtained from BiSON observations (see Chapl in et al 



2002). The solid curve shows the fit consisting of the effect 
of the glitches from the ionization zones and the base of 
the convective envelope, as well as a slowly varying near- 
surface contribution (shown by the dashed line). The lower 
panel shows the individual glitch contributions: the dotted 
and solid lines show the contributions from the first and 
second helium ionization zones, respectively, and the dot- 
dashed line shows the contribution from the base of the con- 



vective envelope. From Hou dek & Gou gh (2007). 



i.e., properties in the star that vary on a scale substan- 
tially smaller than the loca l wavelength of the oscillation s 



dGough & ThompsorJl98d:IVorontsovll98^:lGouphll990b . 
This causes an oscillatory variation in the frequency, de- 
pending on the phase of the eigenfunction at the location 
of the feature; the period of the variation depends on the lo- 
cation of the feature, while its amplitude, as a function of 
frequency, depends on the detailed properties of the feature. 
In the case of solar-like oscillations these effects are related 
to acoustic glitches, with rapid variations in sound speed. 
An important example arises in the second helium ioniza- 
tion zone, as a result of the variation in the adiabatic expo- 
nent Ti = ( d In p/d In p) a( j (the derivative being at constant 
specific entropy), and hence in the sound speed, given by 
c 2 = Tip/p. Glitches also occur at boundaries of convec- 
tive regions where the temperature gradient varies rapidly, 
the sound-speed variation possibly enhanced by variations 
in composition. 



To isolate the effects of the acoustic glitches slower vari- 
ations in the frequencies must be eliminated. An effective 
way of doing that is to consider the second differences, as a 
function of mode order: 

A 2 ^ ti ; = V n -ll - ZVnl + Vn+1 1 (18) 

(e.g.. lGoughlll990|) or possibly higher differences; obvi- 
ously, it is important also to c onsider the error proper- 
ties of the resulting diagnostics. iHoudek & Goughl {2007) 
made a careful analysis of the second difference of 
solar low-degree frequencies, taking into account the 
hydrogen and helium ionization zones and the base 
of the convective envelope; their fit, clearly separating 
these contributions, is shown in Fig. [2] Alternatively, 
one can consider suitable filtering of the ob servations 
( Perez Hernandez & Christensen-Dalsgaardll 1994a ) or sub- 
traction of a fitted smooth function from the freq uencies 
(Mo nteiro et al ] ll994tlMonteiro & Thompsonll2005l) . 

The effects of the second helium ionization obvi- 
ously depend on the helium abundance in the stel- 
lar convection zone. This has been used to deter- 
mine the solar env elope helium abundan c e from helio- 
seismic data (e.g.. IVorontsov et al. 199 lb lAntia & Ba"sul 
1994; Perez Hernandez & Christensen-Dalsgaard 1994b). It 
was sh own bv IPerez Hernandez & Christensen- D al s g aard 
(1998) and iBasu et al.1 (120041) that a similar analysis is in 
principle possible for other main-sequence stars, the ef- 
fect increasing with increasing stellar mass. The proper- 
ties of the base of the solar convective envelope have 
been analysed in considerable detail using the resulting 
acoustic glitch, with emphasis on determ ining the na- 
ture of possible convective overshoot (e.g.. | Basu & Aritial 
1994 iMonteiro et alJll994t i Roxburgh & Vorontsovl 1 1994; 



Christensen-Dalsgaard et al.l 1201 1") . A similar analysis is, 
at least in principle, possible on the basis of stellar 
observations of low-degr ee modes ( Montei ro et al. 2000; 



Mazumdar & Anfiall2001l) . Ballot et alJ (120041) made a de- 
tailed analysis of the potential of determining the depth 
of stellar convection zones from t he frequency vari ations 
caused by the resulting glitch. Also. lMazumd ar (2005) anal- 
ysed the diagnostic potential of the signatures of acoustic 
glitches. 

The effects on the frequencies of acoustic glitches are 
quite subtle and hence it is only with the extensive data from 
CoRoT and Kepler that it has been possible to identify them 
for solar-like oscillations in stars other than the Sun . Base d 
on analysis of CoR oT observations by lCarrier et al. (2010), 
Migli o et al.l (2010) found the signature of the glitch caused 
by the second helium ionization zone in a red giant. They 
used the inferred depth of the ionization zone, together with 
other asteroseismic parameters, to obtain a purely seismic 
determination of the mass and radius of the star; however, 
the data were not yet of sufficient quality to allow a determi 



nation o f the envelope helium abundance. iMazumdar et al. 



d2012al) analysed CoRoT observations of a main-sequence 
star slightly more massive than the Sun and identified a 
very clear signature of the second helium ionization zone 
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and a weak and somewhat ambiguous signature of the base 
of the conve ctive envelope. Very enc ouraging results were 
obtained by iMazumdar et al.1 (12012bl) for several stars ob- 
served with Kepler. 

There may be circumstances where the effects of acous- 
tic glitches might hide or significantly affect other aster- 
oseismic diagnostics. Thus iHoudek & Goughl d201 ll) 'cor- 
rected' solar data for the effect of the glitches associated 
with the second helium ionization zone and the base of 
the convective envelope before using the resulting glitch- 
free frequencies to infer the seismic age and heavy-element 
abundance of the Sun, from an asymptotic analysis of the ef- 
fect of the stellar core. Of course, the by-products obtained 
in such a glitch analysis are also of substantial diagnostic 
value ! 

Sharp features in the sound speed and other seismically 
relevant quantities are produced by the composition struc- 
ture resulting from the evolution of stars with convective 
cores. This is particularly dramatic in moderate-mass stars 
where the mass contained in the convective core grows with 
age, resulting in a composition (and h ence sound-speed) 
discon tinuity if diffusion is neglected. ICunha & Metcalfe 
(120071) derived a diagnostic applicable to this case and 
showed that it provided a measure of stellar age. This 
and similar diagnostic qu antities were further analysed by 
Cunha & Brandao ( 201 lb . A sim ilar analysis was carried 
out bv lsilva Aguirre et al. ( 201 1 ). demonstrating the poten- 
tial of such diagnostics to constrain the size of the convec- 
tive core, including the effects of semiconvection which re- 
main a serious uncertainty in the modelling of stellar evolu- 
tion. 

With the continuing operations of both CoRoT and Ke- 
pler and the resulting longer timeseries there is clearly a 
substantial potential for this type of analysis. The depth of 
convective envelopes is an important parameter in the mod- 
elling of stellar dynamos assumed responsible for stellar 
magnetic cycles, while independent determinations of the 
helium abundance for a range of stars would obviously be 
very valuable for studies of the Galactic chemical evolution 
and as a constraint on Big-Bang nucleosynthesis. 

4 Inversions 

The analysis in the preceding section generally character- 
izes the stellar properties in terms of a limited set of param- 
eters. In particular, in fits such as the minimization of \ 2 (cf. 
Eq.[T7l) the model is characterized by the parameters {Pk}- 
Such fits implicitly assume that the physics of the modelling 
is correct, which is obviously doubtful. Indeed, perhaps the 
most interesting aspect of the asteroseismic analysis is to 
falsify this assumption and, ideally, determine how the mod- 
elling should be improved. 

In the simple minimization of xt an inconsistency in 
the modelling is implied if the minimum of xt is substan- 
tially bigger than unity. However, this is clearly only sig- 
nificant if in fact the errors in the observations have been 



estimated correctly, which in itself is not a trivial task. A 
single mode with an erroneously low estimated error (or a 
single misidentified frequency) could lead to an unrealisti- 
cally large value of xt- Thus care, and perhaps some con- 
servatism, is required in the interpretation of the results of 
the fits. 

If indeed there is convincing evidence that the large 
value of x 2 is significant, the task is obviously to determine 
the cause of the discrepancy, in terms of the stellar mod- 
elling. This would benefit greatly from the ability to locate 
the cause of the frequency differences between the star and 
the best model. Independently of these issues of model fit- 
ting, it is clearly of interest to obtain information about the 
stellar structure that is not tied to specific models, as exem- 
plified by the glitch analysis discussed in Section [3T3l Here I 
provide a brief discussion of inverse analyses with this aim. 

Such techniques have been extremely successful in he- 
lio seismic analysis of solar oscillation data (see, for exam- 
ple |GoughetaL 1 11 9961: Christensen- Dalsgaardll20bl , for re- 
views). Here the mode set is sufficiently rich that it is pos- 
sible to infer the sound speed and density in a very large 
fraction of the solar interior; although this is typically based 
on determining corrections to a reference model, the re- 
sults are generally insensitive to the precise choice of model 
(Bas u et al . 2000). The analysis is based on the assumption 
that the model is sufficiently close the true solar structure 
that the frequency differences 5v n i between the Sun and the 
model can be obtained from an expression linearized in the 
differences in structure^ 



K 



nl 



S r c 2 



5 r p 
P 



dr 



(19) 



Here S r c 2 and S r p are the differences, at fixed r, between 
the Sun and the reference model in squared sound speed 
and density, and Q is a function that accounts for the near- 
surface problems of the model, with a scale factor Q^ 1 that 
depends on the inertia of the modes. The kernels K™\ and 
K™ 1 2 are computed from the eigenfunctions of the modes 
of the reference model. 

The inverse analyses, bas ed on techniques origina lly de- 
veloped in geophysics (e.g. jBackus & Gilbertl 1968), con- 
sist of combining relations of the form given in Eq. (fT9b 
in such a way as to obtain an estimate of, for example, 
6 r c 2 /c 2 that is localized near point, r = rrj, in the Sun. 
If this can be done throughout the star an estimate is ob- 
tained of the sound-speed difference between the Sun and 
the model. Given the extremely accurate solar oscillation 
frequencies spanning degrees from to at least 200 this 
can in fact b e done with high precision between 0.07i? and 
0.95i? (e.g. jBasu et al.lll9971) . The inferences are obtained 
as linear combinations of the frequency differences and can 
be expressed as averages of the true sound-speed difference, 
weighted by averaging kernels /C c 2 „(rrj, r). The success of 



4 Formally, an expression of this form follows from the fact that adia- 
batic oscillations satisfy a variational principle. 
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the inversion is characterized by the properties of PC, par- 
ticularly its width, and the error in the inferred differences. 
The properties of the inversion are typically determined by 
parameters that, for example, control the balance between 
localizing JC and minimizing the errors. The choice of these 
trade-off parameters, and ot her aspects of the helioseism ic 
inversions, was discussed bv lRabello-Soares et al. dl999h . 

In the stellar case we are typically, for the foresee- 
able future, restricted to modes of degree I = 0—3 
and hence the potential for inversion is much more re- 
stri cted. Brief reviews of the related issues we re pro vided 
by Thompson & Christensen-D alsgaard (2002) and Basul 
(120031) . concentrating on solar-like stars. Given the limited 
information it is advantageous to choose a description of 
the stellar interior that, as far as possible, depends on a 
single function of position in the star, while still recogniz- 
ing that the modes observed are predominantly of acous- 
tic nature. Using that c 2 = Y\pj p a reasonable choice is 
the squared isothermal sound speed u = p/p. In fact, in 
much of solar-like stars Ti ~ 5/3 and hence the acous- 
tic properties are fully characterized by u. Assuming fur- 
thermore that the equation of state of the stellar matter is 
known with sufficient accuracy, I?i can be found as a func- 
tion Ti(p, p, {Xi}), where {Xi} is the composition. Char- 
acterizing the dependence on composition simply by the he- 
lium abundance Y by mas^l we can rewrite Eq. ( fT^l l as 



5u n 



K 



111 



S r u 



Kfj r Y 1 dr 



Vnl Jo \ U 

+Q; i lG(Vnl) , (20) 

where in addition we used the equations of hydrostatic equi- 
librium and mass to express the dependence on 8 r p and S r p 
in terms of 8 r u. 

A second important difference between helioseismic 
and asteroseismic inverse analysis is that for stars, unlike 
the case for the Sun, we do not have accurate independent 
determinations of the mass and radius. To separate the in- 
verse analysis from the determination of M and R we can 
rewrite Eq. ( 1201 in terms of a; = r/R, the dimensionless 
function u = {R/GM)u and the dimensionless frequencies 
a nl = 2Tr{R 3 /GM) 1 / 2 is nl : 

SlJnl f /„„7 8 r U 



+ Q; i lG{Onl) ■ (21) 

Here I somewhat crudely assumed that the differences in R 
and M between the star and the reference model are ab- 
sorbed in Q. 

A final difference between the helioseismic and astero- 
seismic case concerns the linear approximation underlying 
Eq. dl9) . It is perhaps reasonable, not least from the tight 
constraints on mass, radius, age and luminosity, to assume 
that solar models are so close to the solar structure that the 
linearized representation of the differences is sufficiently 



accurate. This is not as obviously the case for other stars. In 
fact, analyses by Christensen-Dalsgaard & Thompson (in 
preparation) indicate substantial departures from linearity 
for frequencies of models differing in mass only by a few 
percent from the reference model. These issues need further 
investigation before reliable inversions using this technique 
can be carried out. 

There is relatively limited experience with this type of 
inversion as applied to solar-like stars; in particular it is 
only with the recent space-based data that sufficiently ex- 
tensive data are available to make it realistic to consider 
such analyses. Early tests of inver sion of artificial data 



for low -de gree modes were m ade bv lGough & Kosovichev 



5 In most cases the heavy-element abundance has a negligible effect on 
the thermodynamic state. 



dl993alfbh . Bas u et al.l (12002b tested inversions for model 
frequencies for reasonably realistic sets of stellar artificial 
data, with both models having solar mass and radius. Based 
on Eq. d2~n > they found that relatively localized averaging 
kernels could be constructed in the stellar core, and that the 
inferred S r u showed clear evidence for the core mixing im- 
posed in the test model. An interesting variant of the lin- 
earized inver sion, used to i nfer th e stellar mean density, was 
discussed by iReese et al.l d2012l) who applied it to several 
sets of observed data. 

A conceptually very different technique for asteroseis- 
mic determination of stellar structure has been developed 
by Ro xburgh & Vorontsov (|Roxburgh & Vorontsovl 12002 , 
l2003bh (see also iRoxburghluO 1 01) . This is based on the fact 
that acoustic modes in the outer parts of the star can be char- 
acterized by a phase function which is independent of the 
degree of the mode. By imposing this constraint, at the ob- 
served frequencies, it is possible to determine the structure 
of the bulk of the star as parameterized in terms of a suitable 
variable, such as the density, with fitting parameters that are 
determined iteratively; it is assumed that in the relevant part 
of the star, Ti ~ 5/3. Thus the technique does not depend 
directly on a reference model, although an initial model is 
used as a starting point for the iteration. Applications to arti- 
ficial data have yielded results remarkably close to the orig- 
inal model on which the data were based. 



5 Concluding remarks 

The advances in stellar inferences over the last few years, 
from the observations made by the CoRoT and Kepler mis- 
sions have been remarkable, far beyond the expectations be- 
fore the launch of these missions. In addition to the quality 
of the data, these advances have been based on the rapid 
development of techniques for the analysis of the data. It 
is probably fair to say that the asteroseismic community 
was not entirely prepared for the overwhelming quality, and 
quantity, of the data that we have obtained. 

While striking results have been obtained, we are not 
yet quite at the point of challenging the modelling of stellar 
structure and evolution, with the exception of the evolution 
of stellar rotation where the theoretical basis is still highly 
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uncertain. Such tests and, one might hope, improved under- 
standing of stellar internal physics are of course a key goal 
of asteroseismology. To reach this point the improvements 
in the data quality and quantity resulting from the exten- 
sion of the missions will certainly be very helpful. Equally 
important, however, will be the continued development of 
the techniques used to analyse and interpret the observa- 
tions. This includes a better statistical analysis of the re- 
sults, which will be required to test whether existing mod- 
els deviate in a significant manner from the observations. A 
key tool in this development will be the application of the 
techniques to artificial data; such test s were carried out be- 
fore the launch o f the missions (e.g.. iMonteiro et al.ll2002 : 



Mazumdar 2005) but have since largely been ignored in the 
flush of real data. Similarly, further work is needed to under- 
stand the differences between results of different stellar evo- 
lution and pulsation codes, as started in the ESTA p roject 
in preparation for the CoRoT mission (Monteiro 2008i Fi- 
nally, it is probably time to consider the application of the 
inverse techniques discussed in Section|4] or at least to con- 
tinue developing and testing these techniques on the basis 
of the understanding of the realistic data properties that has 
resulted from the observations. Also, the potential for inver- 
sion localized in the stellar core is undoubtedly improved 
for stars showing m ixed modes, as was demonstrated by 



Lochard et al.l (120041) in the corresponding case of inversion 



for rotation. This deserves further detailed study. 

Beyond the present missions there are proposals for ad- 
ditional missions combining asteroseismology and the study 
of extra-solar planets, such as the Transiting Exoplanet Sur- 
vey Satellite (TESS) under consideration by NASA and 
the PLAn etary Transit s and Oscillations of stars mission 
(PLATO; ICataiall2009h proposed to ESA. If selected, they 
will extend the asteroseismic investigations to more nearby 
stars in very large parts of the sky, although the duration of 
the observations will be shorter than in the case of longest 
timeseries obtained with Kepler. Also, the Stellar Observa 



tions Network group (SONG) network (e.g. Grundahl et al 



201 1) will yield radial-velocity observations of stellar oscil- 



lations of even higher asteroseismic quality than the Kepler 
data, although for a much smaller sample of stars. Thus it 
is crucial to select the astrophysically most interesting tar- 
gets for SONG, ideally based on the understanding of stellar 
properties obtained from the space missions. Here, also, the 
analysis of realistic artificial data is essential. 
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